Controlling conductance statistics of quantum wires by driving ac fields 
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We calculate the entire distribution of the conductance P(G) of a one-dimensional disordered 
system -quantum wire- subject to a time-dependent field. Our calculations are based on Floquet 
theory and a scaling approach to localization. Effects of the applied ac field on the conductance 
statistics can be strong and in some cases dramatic, as in the high-frequency regime where the 
conductance distribution shows a sharp cut-off. In this frequency regime, the conductance is written 
as a product of a frequency-dependent term and a field-independent term, the latter containing the 
information on the disorder in the wire. We thus use the solution of the Mel'nikov equation for 
time-independent transport to calculate P(G) at any degree of disorder. At lower frequencies, it is 
found that the conductance distribution and the correlations of the transmission Floquet modes are 
described by a solution of the Dorokhov-Mello-Pereyra-Kumar equation with an effective number 
of channels. In the regime of strong localization, induced by the disorder or the ac field, P(G) 
is a log-normal distribution. Our theoretical results are verified numerically using a single-band 
Anderson Hamiltonian. 

PACS numbers: 72.10.-d, 72.15.Rn, 73.21.Hb 



I. INTRODUCTION 

Quantum electronic transport in disordered structures 
has been of fundamental and practical interest since semi- 
nal ideas by Anderson and Motti Recently, transport in 
driven systems has received much attention since new 
phenomena have been found and the application of a 
time-dependent driving field opens more possibilities to 
control the transport in an electronic device. For in- 
stance, it has been proposed the possibility of manipulat- 
ing the electronic transport through molecular structures 
by applying an ac field.— ~— 

Disorder in a sample can be an unavoidable and un- 
wanted ingredient in transport experiments, however, it 
can also be seen as an ingredient to be exploited in order 
to manipulate the transport properties of a system. The 
disorder gives a random character to the transport and 
a statistical analysis naturally emerges. In the absence 
of ac fields, the statistics of quantities such as the con- 
ductance has been widely studied. At present, there is a 
good understanding of the effects of disorder on the sta- 
tistical properties of the conductance for one- and quasi- 
one dimensional systems within a non-interacting elec- 
tron model. In fact, the distribution of conductances is 
known, within a scaling approach to localization, at zero 
temperature and infmitcsimally small applied voltage££ 
for one- and quasi-one dimensional systems the evolu- 
tion of the conductance distribution as a function of the 
length of the sample is described by the Mel'nikov and 
Dorokhov-Mello-Pereyra-Kumar (DMPK) equations, re- 
spectively. For finite temperatures and bias voltages as 
well as higher dimensions, some progress has been made 
on the description of the conductance statistics* 7 - 

In contrast to the well studied problem of electronic 



FIG. 1. Schematic view of a one-disordered system of length 
L driven by a frequency-dependent force. The quantum wire 
is connected to electron reservoirs fii and /i2- 



transport in nondriven disordered systems and in spite 
of the possible interest for applications, remarkably lit- 
tle is known about the statistical properties of transport 
quantities when a time-dependent field is applied; effects 
of energy fluctuations on the current statistics in short 
molecular wires has been studied in Ref . Q , while effects 
of ac fields on the localization length and the conduc- 
tance distribution of a ring driven by a time-dependent 
magnetic flux have been studied in Refs. [| and [13], 
respectively. 

In this work, we calculate the complete distribution of 
the conductance of a disordered quantum wire subject to 
a time-periodic driving field. We obtain the distribution 
from high to low frequencies with different regimes of dis- 
order strength and localization. Our theoretical results 
are compared to tight-binding numerical simulations. 



II. 



TIGHT-BINDING MODEL AND 
FLOQUET-GREEN THEORY 



We start by describing briefly the model of our driven 
quantum wire. The conductance is calculated adopting a 
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Floquet scattering approach to the electronic transport 
problem generalizing the Landaucr-Biittikcr formulation 
to driven systems . 11 ' 12 

We might think of an experimental setup where a dis- 
ordered sample, attached by perfect leads to a reservoir 
of electrons on each side, is subject to a laser beam with 
an angle of incidence perpendicular to the wire and a 
polarization angle parallel to the wire, Fig. 1. We as- 
sume that such a system is described by the following 
time-dependent Hamiltonian 



H(t)=H w {t)+H l + H c 



(1) 



The wire Hamiltonian H w (t) is time-periodic: H w (t) = 
H w (t + T), where T = 2tt/uj, uj being the frequency of 
the applied field. For each term in Eq. ((T|) we have: 
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where Ej is the random on-site energy, distributed uni- 
formly in the interval (-W/2,W/2). The form of the 
external time-dependent field assumes a dipolar approx- 
imation for a monochromatic field so the wavelength of 
the laser light must be longer than the length of the wire. 
As the Hamiltonian is time-periodic, the Floquet theorem 
states that there is a set of solutions |</>"' m (i))-Floquct 
states-to the equation 
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where N is the total number of sites of the wire, A is 
the band width of the nondriven system, and V is the 
field amplitude (a factor of 2 is introduced by computa- 
tional convenience) . The Hamiltonian for the leads H\ is 
modeled as 



fc 



(3) 



where c L , R -. k is the creation operator for an electron at 
the left lead (right) with momentum fc and the coupling 
lead- wire Hamiltonian H c is given by 



H c = ^ v ik^ Lk ci + VNkc' m c N + H.c. 



(4) 



The coupling can be described by the spectral density: 
T N{1) (E) = 2TrJ2 k |^ (L ) fc | 2 <5(£ - e k ). The driving field 
might have an effect on the coupling leads- wire, however, 
this situation can be mapped onto a Hamiltonian of the 
form given by Eq. ([1]) by a gauge transformation^ We 
use the wide-band approximation where we assume that 
the coupling between leads and wire is energy indepen- 
dent. This treatment is justified whenever the conduction 
bandwidth of the leads is much larger than all other rel- 
evant energy scales in our problem. We use it as we are 
interested in the effect of the driving in the conductance 
distribution and not of the details of the coupling to the 
leads. In this approximation the information of the wire- 
coupling can be introduced via two self-energies at the 
first and last sites of the wire (Fi and Tn) transforming 
the wire Hamiltonian to a non-Hermitian Hamiltonian. 
Under the above assumptions and within a tight binding 
approximation, the Hamiltonian associated to our quan- 
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where t a,m — e a + mhuj with m an integer number 
and —hut/2 < Re(e Q ) < tko/2. The Green's function 
G(E, t', t") satisfying [IE - H(t')] G(E, t', t")l = U T (t'- 
t"), where I is the identity operator and 5x{t) is a T- 
periodic delta function, can be written in terms of the 
Floquet states \4>f'°(t)} . Thus the Fourier components 
G k (E) of G(E) can be written as 



G {k} (E) = X 
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We define the de-conductance for an ac-driven quantum 
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where I is the current averaged over one period of the 
driving field. Thus G is an experimentally accessible 
quantity which can be written as a sum of the Floquet 
modes g i - k \E F ): 



G= £ gW{E F ), 



where 



g {k) (E F ) 



k=-c 



TiTN\G[ k ^(E)\ 2 being the transmission 



(9) 



(10) 



with T^(E) 

for electrons from the left to the right lead, similarly 
t ni( E ) = T N^i\G i p 1 (E)\ 2 for electrons from the right to 
the left. In the static case, = for systems with 
time- reversal symmetry, however, in the presence of time- 
dependent fields this is not longer truei 2 -. When both 
transmissions are not equal, a current at zero voltage 
can be induced (pumped current). This pumped current 
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FIG. 2. Conductance distributions P[ln(G)] for strength dis- 
order W = 5, field parameters V — 4 and (a) uj = 35 and 

(b) w = 20. L/l = 4.95 and 5.20 for (a) and (b), respectively. 
P(G) for strength disorder W = 2 and field parameters V = 4, 

(c) cj = 35 and (d) cj = 25. L/Z = 0.8 and 0.83 for (c) and (d), 
respectively. A good agreement between theory (solid lines) 
and numerical simulations (histograms) is seen. 



appears as an offset at zero voltage in the I— V charac- 
teristic of our system. The presence of a pumped current 
does not modify our formulas, as we calculate the slope 
of the I—V curve at zero voltage. In the systems we con- 
sidered the pumped current is very much reduced as the 
length of the systems increases and is, in general, negli- 
gible, in contrast to results by Kaiser et. al. in Ref. [8| , 
where they study the current in shorter molecular wires 
(iV = 6). 

We are thus interested in the statistics of the conduc- 
tance G given by Eq. ©. Our theoretical predictions 
are verified numerically by sampling over different disor- 
der realizations, i.e., over different configurations of Ej 
in Eq. ([5]). The numerical simulations are performed 
according to Eqs. ([7]) and where the Fourier com- 
ponents are calculated using the method of matrix 
continued fractions! 9 ' 13 ' 14 We assume perfect lead-wire 
coupling (ri = Ttv = 1) although different coupling 
strengths can be implemented in our numerical simula- 
tions and theoretical framework. The histograms shown 
in all the figures of this work were obtained from 2000 
different disorder realizations. 



III. RESULTS 

A. High-frequency regime 

We start with the simplest case which turns out to be 
the high frequency regime characterized by frequencies 
hui > (A + W)/2. It is also convenient to begin with this 



regime since it illustrates in a simple way the approach 
to our conductance calculations along this work. 

In the absence of disorder and for high frequencies it 
has been shown that the conductance is proportional 
to the square of the Bessel function Jo(2V L j huj) 
A key point is that this frequency-dependent factor 
Jo(2VL/hio) remains under the presence of disorder^ 
Thus, for a disordered wire we write the conductance 
G as the product 



G= Jl(y)g, 



(11) 



where v = 2VL/huj and g is the dimensionless conduc- 
tance of the disordered wire in the absence of the driv- 
ing field. The calculation of the distribution of G is 
now straightforward since the conductance distribution 
for the static case p(g) is given by the solution of the 
Mel'nikov equation, which is given by£ 



p(g) 
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where yo — arccosh(2/g — 1) and s = L/l, 1 being the 
mean free path, s can be seen as a disorder parameter 
which can be obtained from the disorder-average: s = 
(— In g) . In order to calculate P(G) we just need to make 
the change of variable g — > G in Eq. (fl~2j) , using Eq. (fTT|) . 
We can work with the exact expression (|12j) . however, we 
provide an analytical expression, which can be obtained 
by the saddle-point method, for P{G): 



{ ' G 3 / 2 (Jq{v) - G) 1 / 4 
x exp -(Z/£)acosh 2 (j (i/)/VG\ , 



acosh 



(13) 



where Cq is a normalization constant and the ratio L/l 
is obtained from Eq. (fTTj) as 



L/l = 2\nJ {v) - (lnG). 



(14) 



We notice that the L/l is the only parameter in Eq. (fT3")) 
which can be extracted from the numerical experiments 
using Eq. (|14j) ; in this sense, we have a free parame- 
ter theoretical result. In Fig. [5] we compare Eq. (JT3J) 
with the numerical simulations for different values of dis- 
order, frequencies, and field amplitudes. A good agree- 
ment is seen in all cases. We point out the sharp cut-off of 
the distributions, which is determined by the frequency- 
dependent term Jq(v) in Eq. (jlip . This is a manifesta- 
tion of the phenomenon of coherent destruction of tun- 
neling or dynamical localization by the interference of 
multiple paths due to the rapid oscillations of the exter- 
nal field i 15 ' 16 This strong effect on the statistics suggests 
the possibility of controlling the conductance: by tun- 
ing properly the applied field the conductance through 
the disorder sample can be completely switched off at a 
desired conductance value. We finally remark that the 
complete conductance distribution is determined by the 
frequency-dependent term Jo iy) an d the disorder param- 
eter L/l in Eq. (JT3J 
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\a(G) \n(G) 

FIG. 3. Insets: distribution of Sg = |g0 - gl\/(G) for 
L = 10, W = 5, V = 4, uj = 2, and (G) = 6.7 x 1(T 5 (left); and 
L = 10, W = 2, V = 4, w = 2 with (G) = 1.7 x 10" 4 (right). 
Main frames: P[ln(G)] for the same parameters of the insets. 
Solid lines are given by Eq. (|15l) with (InG) = —14.7 and 
a 2 (InG) = 15.5 (left) extracted from the numerical simula- 
tion. (InG) = —11.7 and cr 2 (lnG) = 6 for the right panel. 

B. Intermediate and low- frequency regimes 

Let us consider first the strong localization regime 
(L >> I) where the analysis of the conductance statistics 
is particularly simple. As we shall see, when L becomes 
of the same order of I the calculations of the distribution 
of the conductance are more involved. 

When we increase the disorder strength and/or choose 
the field parameters (amplitude, frequency) such that the 
system is strongly localized, we observe that the Floquet 
modes are strongly correlated; we show two typical ex- 
amples of this behavior in the insets of Fig. 3, where 
we plot the distribution F(Sg) of difference between the 
first and second transmission Floquet modes 5g (we have 
divided this difference between modes by (G) with the 
only purpose of measuring Sg in units of conductance av- 
erage). We observe that F(8g) is highly concentrated at 
small values of Sg. We point out that in the localized 
regime G « 1 and since G is the sum of the Floquet 
modes, Eq. it is expected that each Floquet mode 
is in the localized regime. On the other hand, we re- 
call that in the static case, the conductance distribution 
follows a log-normal distribution. Thus, the distribu- 
tion of the logarithm of the total conductance is given by 
the convolution of log- normal distributions, i.e., each of 
these distributions corresponds to a Floquet mode. The 
resulting distribution is again a log-normal distribution. 
Therefore we write P(lnG) as 

PQnG) = Cg cxp [—(InG — (In G)) 2 /2cr 2 (ln G)] . (15) 

In Fig. 3 (main frames) we show PQnG) for two differ- 
ent disorder strengths. The average and variance of InG 
are extracted from the numerical data and substituted 
into Eq. (fl~5j) . A good agreement between theory and 
numerics is seen. 

Perhaps the most interesting case is when the system is 
not strongly localized by the ac field and/or disorder. We 
go further with our approach to the problem where the 
statistics of the frequency-dependent conductance fluctu- 
ations are derived from the statistical properties of the 



transport problem of the static case. For instance, we 
have found that the Floquet modes g^ show similar cor- 
relations to the conductance channels of the multichan- 
nel (static) case. For simplicity we shall consider the 
case where two Floquet modes g^ and g^ (g*- -1 -* is also 
considered, statistically, it gives the same contribution as 
gM) give the main contribution to Gw^ 

We recall that the DMPK equation is derived for a 
disordered wire whose length L is much larger than its 
width L y (L y « L), i.e., a quasi-one-dimensional sys- 
tem. Within this limit the diffusion of electrons in the 
transverse direction can be neglected, although several 
transverse modes or channels might be open and con- 
tribute to the conductance. For an energy of electrons 
£ and for a given channel n, the longitudinal energy £ n 
is given by £ n = £ — e n , where e„ = mr/L y , n being a 
positive integer. It is important to remark that within 
the DMPK framework only longitudinal diffusion is con- 
sidered and the information of the number of modes is 
given through the energies of the finite number of open 
channels. In our one-dimensional wire driven by an peri- 
odic ac field, the eigenvalues of the Schrbdingcr equation 
are given by e a ' m = e a + mhw, an expression similar 
to the static case; however, m takes unrestricted values 
(positive and negative) whereas the open channels n are 
finite in the static case. With this argumentation we ex- 
pect that the statistical properties of the time-dependent 
problem can be described through the known statistical 
properties of the static case for a quasi-one-dimcnsional 
geometry. As we have mentioned above, we consider the 
case where two Floquet modes give the major contribu- 
tion to the conductance. We shall thus use the DMPK 
results for two open channels. Within the DMPK frame- 
work, the joint distribution of the variables x n related 
to the conductance g n by x n = arccosh\/l/(7„, where n 
labels the channels (n = 1,2), is given by^ 

p(xi,x 2 ) = exp[u(x!,X2) + V(x±) + V(x 2 )], (16) 

where u(xi,X2) = In [| sinh 2 x\ — sinh 2 x 2 \\x\ — x 2 ,] 2 ] 1 / 2 , 
V(x n ) = -3x 2 /2s + ln|a; 2 sinh2x n | 1 / 2 . Equation (TBI) is 
an approximation to the solution of the DMPK equation 
and it is useful from strong to weak localized regimes^. 
As we have discussed above, when the field is applied, 
the Floquet modes follows a statistics given by Eq. 
(HHJ); in order to illustrate this fact we calculate numer- 
ically the distribution of the Floquet modes difference 
d = .g (0) - + g ( ~ 1} ], F(d), and compare to the 

theoretical distribution calculated from Eq. (fl~6"l) . The 
disorder parameter s(— L/l) in Eq. (|16[) is fixed by fit- 
ting the average (d) to its corresponding numerical result. 
In the insets of Fig. 4, F(d) is plotted for two different 
frequency values. As we can see the trends of the numer- 
ical histograms are well described by the theoretical F(d) 
(solid line). These results for F(d) give us confidence in 
using the solution of the DMPK equation, Eq. (fT6|) . Wc 
then calculate G assuming that it can be expressed as the 
product of a frequency-dependent part, a, and a static 
contribution, g: G = ag, i.e., the frequency-dependent 
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FIG. 4. Main frames: conductance distributions for L = 
10, W = 2, V = 1.0, u = 4, a = 0.003 (top) and u = 5, 
a = 0.068 (bottom). The solid line is the analytic result with 
s — 1.5 (top) and s = 1.35 (bottom). The histograms in bro- 
ken lines correspond to the numerical results when the field 
is switched off. Insets: distribution of d = </ ' — \g"' + g^ -1 '] 
for the same parameters of the main frames; theory (solid 
lines) and numerical calculations (histograms) show a good 
agreement. 



term works as a rescaling factor to the static conduc- 
tance, in a similar manner to Jq(v) in the high-frequency 
regime studied previously. This factor is estimated as the 
ratio a = (G) / (<?) , where (<?) is the conductance average 
when the ac field is switched off. Thus, the distribution 
P{G) is given by 

P(G) = (5[G/a- (l/cosh 2 £i + l/cosh 2 x 2 )]), (17) 



where (•) denotes average performed accordingly to Eq. 
([TfJl) . In Fig. 4 (main frames) we compare the theoreti- 
cal distributions (solid line) as given by Eq. (fT7)) and the 
numerical results (histograms) for two different frequen- 
cies: a good agreement is seen in both cases. For dis- 
tinguishing details of the distributions, we plot P(lnG) 
instead of P(G). We also plot P(lnG) for the static case 
(histograms in broken line) in order to remark the strong 
effect of the ac field. Notice the sharp decay of the distri- 
butions, although it is not as strong as the high frequency 
regime. The distributions are also shifted from the origin 
due to renornralization value of the conductance a in Eq. 
(|17p . The value of a depends on some complex interfer- 
ence pattern and as the results in Fig. 4 show, it can 
vary widely with a small variation in the parameters of 
the model. 

C. Summary and conclusions 

The statistical properties of the quantum transport 
in driven disordered systems by ac fields are not well 
known. Here we have studied the effects of an applied 
ac field on the conductance distribution under different 
conditions by varying the disorder strength and field pa- 
rameters. The statistical properties of the frequency- 
dependent problem, in particular the conductance dis- 
tribution, are obtained from known statistical properties 
of the transport in the static case. We have found of 
special interest in the regime of high frequencies where 
the conductance distribution shows a sharp cut-off which 
suggests the possibility of switching applications by tun- 
ing the ac field, although, other aspects such as heat- 
ing or electrons interactions, which we neglect, might be 
relevant to an accurate description of an experimental 
realization. This cut-off is a manifestation of dynami- 
cal localization or coherent destruction of tunneling due 
to multiple interference. Coherent destruction in time- 
periodic potentials has been seen in cold atoms and semi- 
conductor supper latices 3 It should be of interest to ob- 
serve this effect in quantum transport experiments. 

This work was supported by the Programa Ramon y 
Cajal of MICINN (Spain), the European Social Fund, 
and MICINN (Spain) under Project No. FIS2009-07277. 
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